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models  the  height  change  at  the  ridge  as  an  effective  shunt  susceptance.  Calculated  and  measured 
data  are  compared;  they  show  good  agreement.  The  analysis  should  prove  useful  in  designing 
transformers  to  match  empty  ridged  waveguide  to  slab  loaded  rectangular  waveguide  components 
such  as  toroidal  phase  shifters. 
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TE-MODE  SOLUTIONS  FOR  DIELECTRIC-SLAB 
CENTER-LOADED  RIDGED  WAVEGUIDE 


INTRODUCTION 

Applications  exist  which  require  dielectric  or  ferrite  slab  center  loaded  rectangular 
waveguide  to  be  used  in  conjuction  with  ridged  waveguide.  Transitions  can  be  made  with 
stepped  matching  transformers;  these  transformers  are  appropriate  sections  of  a composite  of 
the  two  different  waveguide  configurations.  The  objective  of  this  report  is  to  present  an  analysis 
of  such  dielectric  slab  center  loaded  ridged  waveguide  and  to  provide  a method  of  obtaining 
equations  for  the  TEll0  propagation  characteristics,  thus  facilitating  transformer  design. 

BACKGROUND 

Because  higher  order  modes  can  easily  cause  mismatch  and  transmission  loss  spikes, 
waveguide  operation  is  generally  limited  to  a frequency  band  where  only  the  principal  mode 
may  propagate.  Conventional  rectangular  waveguide  has  a theoretical  two-to-one,  or  sing’e  oc- 
tave, principal-mode-only  frequency  bandwidth;  in  practice,  the  useable  bandwidth  is  lesj  be- 
cause of  large  attenuation  near  the  cut-off  frequency. 

Ridged  waveguide,  particularly  double-ridged  waveguide,  is  commonly  used  when  larger 
bandwidths  are  required  at  high  power  levels.  A frequency  range  of  more  than  four  to  one 
between  the  cut-off  frequencies  of  the  TEl0  and  TE20  modes  can  easily  be  obtained  [1,2]  using 
double-ridged  waveguide.  Similar  bandwidths  can  be  obtained  with  dielectric  slab  center  loaded 
rectangular  waveguide  [3,4,5].  Both  ridged  and  slab  loaded  waveguide  achieve  broad 
bandwidths  by  adding  large  capacitance  to  the  dominant  mode  while  only  slightly  affecting  the 
capacitance  of  the  next  higher  order  mode. 

Ferrite  toroidal  phase  shifters  also  can  be  designed  for  operation  in  excess  of  one  octave. 
Because  of  the  small  gap  spacing  of  ridged  waveguide,  the  phase  shifters  are  generally  made  in 
rectangular  waveguide.  Dielectric  slab  center  loaded  rectangular  waveguide  would  be  readily 
compatible  with  the  ferrite  toroidal  phase  shifter,  but  it  is  not  a commonly  used  transmission 
line.  Since  ridged  waveguide  is  commonly  used,  it  would  be  desirable  to  have  compatibility,  i.e., 
matching  transitions,  between  ridged  waveguide  and  ferrite  toroidal  loaded  rectangular 
waveguide.  Dielectric  loaded  tapered  transitions  are  possible,  but  the  fabrication  would  be  very 
difficult.  Also,  a quasi-Tchebycheff  transformer  design  should  give  better  matching  for  given 
length  transitions.  The  latter  approach  requires  transformer  sections  of  dielectric  loaded  ridged 
waveguide,  but  analysis  of  this  type  of  transmission  line  is  not  currently  available  in  the  litera- 
ture. The  analysis  in  this  report  employs  an  equivalent  transmission  line  circuit  for  the 
transverse  component  of  the  propagating  electromagnetic  wave  to  derive  solutions  for  the  TEllQ 
propagation  constants  in  such  waveguide. 
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ANALYSIS 

TE-mode  solutions  for  the  dielectric  slab  center  loaded  rectangular  waveguide  of  Fig.  la 
can  be  derived  by  using  ABCD  matrices  [6]  or  by  using  an  equivalent  transmission  line  circuit 
for  the  crossguide  component  of  the  electromagnetic  wave.  The  homogeneous  double-ridged 
waveguide  of  Fig.  lb  has  been  analyzed  [1,2]  by  using  the  latter  method  in  conjunction  with 
the  equivalent  discontinuity  susceptance  due  to  the  height  change  at  the  ridge  wall. 


(a)  Dielectric  slab  center  loaded  rec-  (b)  Double-ridged  waveguide 

langular  waveguide 


Fig.  I — Broadband  Waveguide  cross  sections 

For  the  dielectric  slab  center  loaded  double  ridged  waveguide  of  Fig.  2,  the  analysis  is 
similar  to  that  for  the  homogeneous  case,  with  an  extra  section  incorporated  in  the  equivalent 
transmission  line  circuit.  The  dimensions  referred  to  in  all  subsequent  calculations  are  those 
shown  in  Fig.  2.  For  simplicity,  this  report  will  consider  only  the  case  for  TEn0  modes  and  will 
assume  that  the  transmission  line  is  lossless,  i.e.,  perfectly  conducting  waveguide  walls  and  a 
dielectric  loss  tangent  of  zero.  Axial  symmetry  will  also  be  assumed. 


Fig.  2 — Cross  section  of  dielectric  slab  center  loaded, 
double-ridged  waveguide 


Cohn’s  article  on  ridged  waveguide  [2]  points  out  that  for  the  homogeneous  case  (i.e.  «f  = 
1)  the  cross  section  may  be  treated  at  the  cut-ofT frequency  by  assuming  that  it  is  an  infinitely 
wide,  composite,  parallel  strip  transmission  line  short-circ  .ted  at  two  points.  The  resultant 
electromagnetic  field  may  be  considered  as  an  electromagnetic  wave  traveling  from  side  to  side 
without  longitudinal  propagation.  The  resonant  conditions  can  then  be  solved  for  the  cut-off 
frequencies  of  the  different  TEll0  modes. 

A similar  argument  holds  for  the  inhomogeneous  case.  In  addition,  the  longitudinal  pro- 
pagation constant  may  be  treated  as  the  unknown  quantity,  and  solutions  at  any  frequency  may 
be  obtained  by  separating  the  wave  vector  in  each  region  into  its  transverse  and  longitudinal 
components.  Since  the  waveguide  configuration  is  symmetrical,  the  resonance  condition  for  the 
transverse  wave  component  will  result  in  an  infinite  (zero)  impedance  at  the  center  for  n odd 
(even).  Half  of  a cross  section  is  shown  in  Fig.  3a,  and  the  equivalent  transmission  line  circuits 
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for  the  transverse  wave  are  shown  in  Fig.  3b  for  n odd  and  Fig.  3c  for  n even.  Since  the 
equivalent  circuit  is  a composite,  dissipation  less,  passive  line  matched  at  both  ends,  it  is 
matched  at  all  points.  Therefore,  the  sum  of  the  admittances  at  the  plane  y2  of  the  effective 
lumped  capacitance  due  to  the  ridge  wall  must  equal  zero.  Within  each  region,  where  the  re- 
gions are  shown  in  Fig.  3a,  2^,  is  the  characteristic  impedance,  >o,  = 1 /Zq,  is  the  characteristic 
admittance,  and  0,  is  the  transverse  electrical  length;  0,  is  equal  to  the  product  of  the  physical 
transverse  dimension  of  the  region  and  yyn  the  complex  transverse  propagation  constant.  Since 
all  regions  are  lossless,  y , and  therefore  0,,  will  be  purely  real  or  purely  imaginary. 

region  (1)  (2)  (3),  Electric  (n  even) 

” or  magnetic  (n  odd) 

T . I wall 
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(b)  Equivalent  transmission  line  circuit 
for  » odd 


(c)  Equivalent  transmission  line  circuit 
for  n even 


Fig.  3 Half  waveguide  cross  section  and  equivalent  transmission  line  circuits 
for  transverse  wave 


The  reflected  impedance  Z presented  by  a load  impedance  ZL  terminating  a transmission 
line  of  characteristic  impedance  Zq  with  propagation  constant  y and  length  w is  [7] 

(zL  + z0)ey- + (zL  - z0)e~y- 

Z = Z V I ) 

^ (ZL  + 2^)eyw  - ( ZL  -z^)e-yw 
The  short  circuit  at  in  Fig.  3b  will  be  reflected  back  to  y2  as  ^1-2  where 


zl-2  “Zoi,anh  y^i  j 


^1-2  “ >01  co,h  I K 


a — s 
2 
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The  open  circuit  at  y4  will  reflect  back  to  y3  as  Z4_3  with 


Z4-3  “ 3)3  coth  1^3  yj 

Equation  (1)  can  be  rewritten  in  the  form 

ZL  cosh  yw  + Zq  sinh  yw 
ZL  sinh  yw  + Zq  cosh  yw' 

Since  Z4_3  terminates  region  2, 

Z4_3  cosh  yy2  S 2 ' + Zq2  sinh  yy2  5 ^ ' 

Z4-2  = 3)2  ■■ 

Z4_3  sinh  yy2  S-  + Zq2  sinh  yy2  S 

Using  6 i = yyiw,  to  simplify  notation  and  substituting  Eq.  (4)  into  Eq.  (6)  yields 
£ Z^j3  coth  03  cosh  02  + Zqj  sinh  02 

Z4-2  3)2  ^ coth  g ^ sjnh  2^2  cosh  02 


' Zqj  coth  Oj  sinh  02  + Zq2  cosh  02 

4-2  02  Zqj  coth  Oj  cosh  02  + Zq j sinh  02 

Since  the  sum  of  admittances  at  y2  must  equal  zero, 

3-2  +JBC  + K4_2  -0. 

Substituting  Eqs.  (3)  and  (8)  into  Eq.  (9)  yields 

Zqj  coth  Oj  sinh  0,  + Zq2  cosh  0, 

Yq , coth  0,  +JB.  + Yq2  22 — f — J1 2.  2 

Ul  I r 02  ^ COth  g 3 cosh  g 2 + Zq 2 sinh  02 


M)2 

coth  03  sinh  02  + — — cosh  02 
Bc  '02  M)1 

coth  0,  + y 22 =0.  (11) 

*01  '01  ^-02 

coth  03  cosh  02  + sinh  02 

3)3 

Since  region  1 and  region  2 have  the  same  propagation  constant,  yyi  = yy2,  the  impedances 
are  proportional  to  the  heights: 

(12) 

3).  »02  * 

Regions  2 and  3 have  equal  heights,  and  since  the  transverse  wave  is  TE , the  impedance  ratio 
is 

£--2*-.  (13) 

3)3  Tji2 

The  left  side  of  Eq.  (11)  may  be  rewritten  as  a single  fraction.  All  terms  in  the  denominator  are 
finite,  so  the  numerator  may  be  equated  to  zero.  The  resultant  expression  is 
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fl 


— sinh  9, 
a 


lyv2  cosh  0 3 sinh  02  + yy3  s'0*1  e3  cosh  1 + 


cosh  0|  + j ■ 


B. 


'01 


sinh  0 


x [-yi2  cosh  03  cosh  9 2 + yy$  sinh  03  sinh  02  ] = 0. 


Within  each  region. 


where 


'I  + yji  + yl  = _<u  W 


= «0  for  / =1,2 

tre0 


(14) 


(15) 


For  TE  modes,  yxi  = 0 for  all  regions  and  y,,  = ji 3 for  all  regions;  /3  is  the  longitudinal  propa- 
gation constant  (above  cutoff)  for  the  waveguide  configuration.  Substituting 

, n 2 


yyi  = y/p2 

2 M0  «,  for 

Q> 

= j-J(o  2 Mo  « 

i - P2  for 

to 

and 

= yy,'v, 

with 

"l 

= 1/2  (a  - 

s) 

w2 

= 1/2  (s  - 

t) 

*3 

= 1/2  / 

2 e,  > P2 


(16) 


(17) 


into  Eq.  (14)  yields  the  transcendental  equation  in  p that  must  be  solved  for  TEn0  (n  odd) 
modes.  The  smallest  root  of  Eq.  (14)  is  the  TEi0  solution,  the  next  root  the  TEi0  solution,  etc. 

For  TEn o (weven)  modes,  the  analysis  starts  with  the  equivalent  transmission  line  circuit 
of  Fig.  3b  and  proceeds  in  a manner  similar  to  the  case  for  n odd.  The  resultant  transcendental 
equation  is 

' B, 

[■y,2  sinh  03  sinh  02  + yvj  cosh  03  cosh  02  ] + 


— sinh  0i 
d 


cosh  0 , + j ■ 


sinh  0 


'01 


x [yv2  sinh  03  cosh  02  + yyi  cosh  03  sinh  02  ] =0 
with  Eqs.  (16)  and  (17)  being  applicable. 


(18) 


If  tr  = 1,  it  is  straightforward  to  show  that  Eqs.  (14)  and  (18)  reduce  to  the  expressions 
for  the  odd  and  even  mode  cutoff  frequencies,  respectively,  for  double-ridged  waveguide  11,2]. 
Also,  if  b = d,  Bc  equals  zero  and  Eqs.  (14)  and  (18)  result  in  expressions  for  the  odd-  and 
even-mode  propagation  constants  of  dielectric  slab  center  loaded  rectangular  waveguide  identi- 
cal to  those  obtained  by  use  of  ABCD  matrices  16]. 

The  discontinuity-susceptance  term  S(/T01  is  obtained  from  the  Waveguide  Handbook 
[81.  Appendix  A gives  the  necessary  equations  for  calculating  Bc/Y0 ( in  terms  of  the 
waveguide  dimensions  (from  Fig.  2)  and  the  effective  wavelength  A^,.  Note  that  AR  is  the 
wavelength  of  the  wave  component  which  is  incident  normal  to  the  height  change.  Therefore 
A of  Appendix  A is  the  wavelength  of  the  transverse  wave  in  regions  I and  2,  namely  Av|. 
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For  standard  (i.e.  air  filled)  double-ridged  waveguide  er  =1;  thus 

?yl  = yy2  =>>-3  =JPy  V 

and  X^,,  = 2n/fiy  is  a real  constant  for  a given  configuration.  However,  for  the  general  case 
tr  > 1 and 

y>  \ = y v2  * y,  3 ■ 

with  the  result  that  the  values  of  yyi  that  satisfy  Eq.  (14)  or  (18),  subject  to  (16)  and  (17),  are 
no  longer  constant  but  depend  on  the  frequency.  This  is  of  course  to  be  expected,  since  for 
any  nonhomogeneous  waveguide  a l/vl  — (fc//) 2 term  no  longer  describes  the  dispersive  na- 
ture. However,  there  is  another  problem  because  of  the  inhomogeneity.  At  all  frequencies 
above  cutoff-,  the  transverse  wave  propagation  constant  in  the  dielectric  region  will  be  entirely 
imaginary;  i.e., 

yyi  = JPyi  for  " > 

However,  there  is  a critical  frequency  («»,.„,  is  greater  than  the  cutoff  frequency  how 
much  greater  depends  upon  the  degree  of  dielectric  loading)  such  that  for  frequencies  greater 
than  atcrll  the  transverse  propagation  constant  in  regions  1 and  2 is  real,  that  is, 

yy l = Yyl  ~ a/l  for  > ",  nr 

When  (i)  > <ocnl,  the  transverse  "wave”  in  these  regions  is  no  longer  a resonant  traveling  wau 
but  rather  the  fields  are  decaying  exponentially  away  from  the  dielectric  region,  and  the  con- 
cept of  wavelength  in  the  region  of  the  discontinuity  is  not  meaningful.  The  expression  for  the 
Bc/y0i  term  from  Ref.  8 is  no  longer  applicable;  indeed,  the  validity  of  the  equivalent 
transmission  line  circuit  for  the  waveguide  height  change  (a  shunt  susceptance  at  the  junction 
of  two  transmission  lines  of  unequal  characteristic  impedance)  is  questionable  for  operation 
below  cutoff.  Also,  the  calculation  for  Bc/Y0i  is  based  on  a model  which  assumes  that  the 
waveguide  extends  to  infinity  in  both  directions  away  from  the  height  discontinuity;  in  prac- 
tice, the  assumption  is  valid  if  additional  mismatches  are  far  enough  removed  from  the  height 
discontinuity  so  that  the  local  fields  have  decayed  to  small  proportions.  These  local  fields  are 
the  evanescent  modes  of  the  fringing  fields  caused  by  the  height  discontinuity,  and  they  decay 
very  rapidly. 

Future  investigation  is  planned  to  model  an  equivelent  circuit  of  the  waveguide  h#ght 
change  to  include  operation  below  as  well  as  above  the  cut-off  frequency,  and  to  include  the 
proximity  effects  of  waveguide  walls  and  dielectric  center  loading.  However,  for  this  report  the 
following  two  engineering  assumptions  are  made:  9 


I.  The  Bt./Y0i  term  can  be  neglected  for  frequencies  below  the  critical  frequency.  Si:  ce 

B, 


0 as  <d  — > * 


'ot 


and  for  <o  < wrril  the  fields  of  the  transverse  wave  are  decaying  exponentially  in  the  region  of 
the  height  discontinuity,  a small  shunt  susceptance  term  will  have  only  a minor  effect  on  the 
solution  for  fi.  Equations  (14)  and  (18)  are  transcendental  equations  and  must  be  solved  by 
some  algorithm  using  trial  values  of  fi-  If  a trial  value  of  fi  yields  an  imaginary  transverse  pro- 
pagation constant  in  region  1,  the  BJ  *01  term  is  calculated  with 

Yy\  ~ JPy\  and  X vj  =— . 

P y\ 
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If  the  trial  value  of  /3  results  in  real,  the  Bc/  n>i  term  is  neglected,  i.e.  set  equal  to  zero,  in 
Eqs.  (14)  and  (18). 

2.  Proximity  effects  can  be  neglected  in  the  calculation  of  BJ  *01  • 

Although  the  validity  of  these  two  assumptions  may  be  questioned  from  a rigorous 
theoretical  aspect,  the  close  agreement  between  calculated  and  measured  values  of  fi  for 
different  configurations  (shown  in  Figs.  4 and  5)  indicates  that  both  assumptions  result  in  accu- 
racy sufficient  for  most  practical  applications. 

A listing  of  a computer  program  to  solve  for  the  principal  (TEW  ) mode  propagation  con- 
stant of  dielectric  slab  center  loaded  double-ridged  waveguide  is  given  in  Appendix  B. 

All  discussions  and  calculations  thus  far  have  assumed  a double-ridged  waveguide 
configuration.  For  the  asymmetric  or  single-ridged  waveguide  configurations  shown  in  Fig.  6, 
Eqs.  (14)  and  (18)  remain  valid;  however,  the  expression  for  Bc/  *01  must  have  Xv,[  replaced 
by  1/2  X^, . 


^(degrees/inch)  ({(degrees/cm) 


75  8 9 TO  ff  12  13  t4  ?5  ?6  17  18 


Frequency  (GHz) 

Fig  4 — Calculated  and  measured  values  of  fi  for  dielec- 
tric slab  loaded  double-ridged  waveguide 
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pMegrees/mch) 


..((ivyrees/cm) 


Fig.  5 — Calculated  and  measured  values  of  /3  for  dielec- 
tric slab  loaded  doublcd-ridged  waveguide 


Fig.  6 — Cross  section  of  dielectic  slab 
center  loaded  single-ridged  waveguide 


CONCLUSION 

Based  on  the  equivalent  transmission  line  circuit  for  the  transverse  component  of  the 
propagating  electromagnetic  wave,  expressions  have  been  derived  for  the  TEn 0 mode  propaga- 
tion constants  of  a dielectric  slab  center  loaded  ridged  waveguide  configuration.  These  expres- 
sions are  transcendental  equations  involving  the  propagation  constant,  but  they  can  readily  be 
solved  with  a computer.  Based  on  the  agreement  between  calculated  and  measured  data,  cer- 
tain assumptions  made  in  the  derivation  appear  valid.  The  analysis  should  prove  useful  in 
designing  transformers  to  match  ridged  waveguide  to  dielectric  or  ferrite  slab  center  loaded  rec- 
tangular waveguide. 
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Appendix  B 

FORTRAN  LISTING  OF  COMPUTER  PROGRAM 


00100 

C THIS 

IS  PROGRAM  DRWGDL. FOR  - CWY  -OCT  75 

00200 

INTEGER  RIK 

00300 

REAL  k'XAIR 

00400 

PI=3. 1415927 

00500 

C=2. 997925E+08 

00600 

Rl=39. 37003 

00700 

R2=2. 0»R1 

00300 

RRMDI = 180. 0/<PI*Rl) 

00900 

Cl  = (2. 0E+09*PI/C)**2 

01000 

NEWRUN=  0 

01100 

TYPE  600 

01200 

600 

FOPMAT  <///'  PROGRAM  DRWGDL  CWY  OCT  75'/"'  COMPUTES 

01300 

1 

TE10  CUTOFF  FREQUENCIES  AND  PROPAGATION  CONSTANTS  OF'/ 

01400 

2 

' SYMMETRIC  DIELECTRIC  LOADED  DOUBLE  RIDGED  WAVEGUIDE') 

01500 

1 05 

TYPE  605 

01600 

6 05 

FORMAT <//'  WAVEGUIDE  DIMENSIONS  IN  INCHES  - A»B*D*S»  ' 

01700 

READ  <5*  ♦)  A j Bj  D j S 

01300 

1 06 

TYPE  615 

01900 

615 

FORMAT  </'  RELATIVE  DIELECTRIC  CONSTANT  OF  CENTER 

02000 

1 

LOADING  = •'$) 

02100 

READ  C5» ♦ > EPSR 

02200 

TYPE  625 

02300 

625 

FORMAT  <'  WIDTH  IN  INCHES  OF  CENTER  LOADING  = '*) 

02400 

ACCEPT  630- T 

02500 

630 

FORMAT  <F3.3‘ 

02600 

IF  <T . LT . S) GO  TD  103 

02700 

TYPE  631 

02300 

631 

FORMAT  <'  DIELECTRIC  WIDTH  MUST  BE  LESS  THAN 

02900 

1 

RIDGE  WIDTH  TRY  AGAIN) 

03000 

GO  TO  105 

03100 

108 

TYPE  606 

03200 

• 606 

FORMAT  <//'  DRWGDL  PARAMETERS  DIMENSIONS  IN 

03300 

1 

I NCHES ' /8X ' A'9X ' B ' 9X  D ' 9X ' S ' 1 2X  T 6X4HEPS '/IX*) 

03400 

TYPE  607.A>B.D»S»T> EPSR 

03500 

607 

FOPMAT  C4F10.4.F13.4.F10.3) 

03600 

R=D/B 

03700 

RS=R*»2 

03300 

IFP=0 

03900 

IF  <ABS  'R-l . 0)  . LT.  1.  0E- 06)  IFP=1 

04000 

Wl= (A-S) /R2 

0*»100 

W2=  (S-T) /R2 

04200 

W3=T'R2 

04300 

CEREST= 1 . + <1 . / R- 1 . ) ♦COS  <P I ♦ ' A-S) / (2 . ♦A) ) 

04400 

CLREST=CEREST+ (EPSR-1 . ) /R*COS  <PI»(A-T) / <2. *A) ) 

04500 

EDCTPY=CLPEST/CEREST 

04600 

ALCEST=A*CLREST» <R+  <1 . -R) ♦SIN  (PI*  <A-S) / <2. *A>  ) ) 

04700 

C THF 

ABOVE  FOUR  QUANTITIES  ARE  TO  BE  USED  FOR  CALCULATING 

04300 

C APPROXIMATE  CSTAPTING  VALUES)  OF  CUTOFF  FREQUENCIES  AND 

04900 

C PPOPAGAT ION  CONSTANTS 

05000 

IBC=1 

12 
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05100 

FREQ=C^R 1 s CALCEST *2 . 0E+09) 

05200 

DELBY=0.31*FREQ 

05300 

BY=0. 0 

05400 

60  TO  112 

05500 

109 

IF  CNEWRUN. LT. 2>  GO  TO  210 

05600 

IF  CFSTART. GT. FCGHZ)  60  TO  220 

05700 

210 

TYPE  635 

05300 

635 

FORMAT  <✓'  FREQUENCIES  IN  GHZ  - START j STOP? INCREMENT : 

05900 

READ (5. ♦>  FSTART  > FSTOP » DELF 

06000 

640 

FORMAT  CF9. 3> 1X»  F9. 3> 1X>F9. 3> 

06100 

220 

IFCFSTAPT.LT. 1. OE-13)GO  TO  180 

06200 

IF  CFSTART. GT. FCGHZ?  GO  TO  230 

06300 

TYPE  645 

06400 

645 

FORMAT  C'  FREQUENCY  MUST  BE  GREATER  THAN  CUTOFF  '> 

06500 

GO  TO  210 

06600 

IFCFSTOP.LT. 1. OE-13>FSTOP=FSTAPT-1. 0 

06700 

230 

TYPE  655 

06800 

655 

FORMAT  C'4X4HFREQ8X4HBET A9X3HGWL7X5HRAT I 08X5HKXA I R/ 

06900 

1 

5X3HGHZ6X6HDEG^IN6X6HINCHES4X8HGWL/FSWL7X6HR  OR 

07000 

110 

IFREQ=0 

07100 

FREQ=FSTART 

07200 

1 1 1 

IFREQ=IFREQ+1 

0?300 

BY=P I *2 . E+ 09/C^SQRT  CEDCTRY* CFPEQ~2+FCGHZ**2>  > 

07400 

C THIS 

IS  A FIRST  TRY  FOP  BETA 

07500 

DELBY=-0.31*BY 

07600 

112 

ICRDSS=0 

07700 

ITAN=0 

07300 

IBTRY=0 

07900 

115 

C1F=C1»FREQ**2 

03000 

C1FEP=C1F*EPSR 

03100 

120 

IBTRY=IBTRY+1 

03200 

IFCIBTRY.LT. 265GO  TO  122 

03300 

TYPE  705 

03400 

705 

FORMAT  MORE  THAN  25  TRIES  AT  ROOT') 

03500 

GO  TO  170 

03600 

122 

BYSQ=BY**2 

03700 

GX3SQ=C 1 FEP-BYSQ 

03300 

GX2SQ=C 1 F-BYSQ 

08900 

6X3=SQPT CABS  CGX3SQ) > 

09000 

GX2=SQRT  CABS  CGX2SQ) > 

09100 

IF  CGX3SQ) 1 3 0 » 132» 132 

09200 

130 

CHS3=SINHCGX3*W3> 

09300 

CHC3=C0SH  CGX3*W3) 

094  00 

IPGX3=1 

09500 

GO  TO  134 

09600 

132 

CHS3=SINCGX3*W3) 

09700 

CHC3=CDS C6X3*M3) 

09300 

IRGX3=-1 

09900 

134 

CONTINUE 

10000 

IF  CGX2SQ) 136» 138> 138 
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10100 
1 0200 
1 0300 
10400 
10500 
10600 
10700 
10300 
10900 
11000 
11100 
11200 
11300 
11400 
11500 
11600 
11700 
1 1800 
11900 
12000 
12100 
12200 
12300 
12400 
12500 
12600 
12700 
12800 
12900 
1 3000 
13100 
13200 
13300 
13400 
13500 
13600 
13700 
133  00 
13900 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
15000 


136  CHS2=SINHCGX£*W2> 

CHC2=C0SH  <GX2+W£> 

CHS1=SINHCGX£*W1 > 

CHC 1 =COSH  CGX2+W 1 > 

KXfiIP=GX2 
RIK=1HP 
IRGX2=1 
GD  TO  140 

138  CHS2=SINC6X2*W£> 

CHC2=C0S  CGX£+W2> 

CHS 1 ®S 1 N CGX£»W 1 > 

CHC 1 =CDS  CGX2+W 1 > 

IRGX2=-1 
KXfi  I R=GX£+PRMB I 
RIK=1HI 
140  BDY=0. 0 

IF  CIFR.EQ. 1>G0  TO  153 
IF CIRGX2. EQ. 1> GO  TO  153 
C CALCULATE  B 'Y  TERM 

p=<i+r:>^<i-r> 

GL=2.  0+PI.'GX£ 

P£=SGRT  C 1 . 0-  CB <R  1 +GL>  > ♦♦£> 

P3=SQRT C 1 . 0- CIV  CR 1 ♦GL> > ♦♦£> 

PA=P+*  <g.  04R)  ♦ *1 . 0+P£>  ✓ Cl . 0-P2>  - c 1 . 0+3.  0+RS> /(l.  0-PS> 
PhF’=P++  C£.  0'R>  * Cl . 0+P3>  s<l.  0-P3>  + <3.  0+RS>  ✓ Cl . 0-PS> 

PC=  C <4 . 0»R)  ■'  ( 1 . 0-RS)  > ♦♦£ 

PT1=AL0G  C <1 . 0-RS>y<4.  0»R> +P++ C 0. 5* CR+1 . 0 R> > > 

PT£*2.  0*  CFA+PAP+2.  0+PO  /•  CPA+PAP-PC++2' 

PT3=  CB  CP  1 +4 . 0+GL>  > ++2+  C 1 . 0-P>  ♦♦  C4 . 0+P>  ♦ C <5 . 0+PS 
1 - 1 . 0> ✓ C 1 . 0-RS> +4 . 0+RS+PC <3 . 0+PA > > ♦♦£ 

B0Y=2. 0+B+ CPT 1 +PT2+PT3> / CR1+GL- 
C CALCULATE  FCBETA.’ 

1 53  FBETA=R+  C-BOY+CHS 1 +CHC 1 ) ♦ C 6 X 2 • G X 3 * C H C 3 ♦ C H C 2 

1 + 1 RGX3+CHS3+CHSS) + 1 PGX2+CHS 1 ♦ CGX2 ■ GX3+CHC3+CH  52 
£ + I RGX2+ I PGX3+CHS3+CHC2 > 

IFCIBC.EQ. 1>  BY=FREQ 
C ROOT  SEARCH  ROUTINE 

IF  CABS  CFBETA1 . LT. 1 . 0E-03> GD  TO  170 
IFCIBTRY.EO. l>GO  TQ  163 
IF  CITAN.EQ. 1>GD  TQ  164 
IFCFBETA+FBOLD.LT. 0. O’GO  TQ  161 
IF  CABS  CFBETfi’ . GT. ABS  CFBQLD’ > DELBY=-BELBY 
GQ  TO  162 

161  DELBY=-DELBY 
ICRQSS=1 

162  IF  CICRQSS. EQ. 1> DELBY=0. 5+DELBY 

163  BYNEUI=BY+DELBY 

IF  CABS C CBY-BYMEM>  , BY> . LT. 0. 1 > I TfiN= 1 
130  TO  166 

164  IF CfiBS CBY-BYQLD> .LT. l.E-05.flNH. FBETH.LT. l.E-06> GO  TQ  170 
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15100 

BYNEW=BY-FBETA»  <BY-BY0LD> / <FBETA-FBOLD> 

15200 

166 

BYDLD=BY 

15300 

BY=BYNEW 

15400 

FBOLD=FBETfi 

15500 

IF  (IBC. EQ. 2>  GD  TO  120 

15600 

FREQ=BY 

15700 

BY=0. 0 

15800 

GO  TO  115 

15900 

170 

IF  <IBC. EQ. 2>  GO  TO  175 

16000 

FCGH2=BY 

16100 

TYPE  658*  FCGHZ*  BOY 

16200 

658 

FORMAT  </"'  TE1 0 MODE  CUTOFF  FREQUENCY  IN  GHZ  = ' F7.4' 

16300 

1 B/Y  = 'F7.3> 

16400 

IBC=2 

16500 

GO  TO  109 

16600 

173 

CONTINUE 

16700 

175 

BYDI=BY*RRMDI 

16800 

GWL=360. 0/BYDI 

16900 

FSWL=R1*C/<FREQM. 0E+09’ 

17000 

RGLFS=GWL^FSWL 

17100 

177 

TYPE  66 0 * FREQ*  BYB I * 6WL * RGLFS * K.XR I R * R I K 

17200 

660 

FORMRT  <1X*F7.3*3X*F9.2*3X*F9.4*4X*F8.4* 3X*F8.2* 1X*R1> 

17300 

IF  <FREQ.  GE.  FSTOP5  GO  TO  130 

17400 

FREO=FREQ+DELF 

17500 

GO  TO  111 

17600 

180 

TYPE  665 

17700 

665 

FORMRT  <✓/"  WISH  NEW  PARAMETERS?  NONE=0*  RLL=1 * 

17300 

1 CENTER  L0ABING=2*  FREQ=3 

17900 

ACCEPT  670*  NEWRUN 

18000 

670 

FORMRT  <I1> 

13100 

GO  TO  099*  105*  106*210*  1805  NEWRUN+1 

13200 

199 

CONTINUE 

13300 

END 

15 


